Statistical Model of Superconductivity in a 2D Binary 
Boson-Fermion Mixture 



M. Casas a , N.J. Davidson 6 , M. de Llano c , T.A. Mamedov d 
A. Puente a , R.M. Quick 6 , A. Rigo a and M.A. Soli's 6 

February 1, 2008 

a Departament de Fisica, Universitat de les Illes Balears, 07071 Palma de Mallorca, Spain 
^Department of Physics, University of Pretoria, 0002 Pretoria, South Africa 
c Instituto de Investigaciones en Materiales, Universidad Nacional Autonoma de Mexico, 
Apdo. Postal 70-360, 04510 Mexico, DF, Mexico 

^Department of Physics Engineering, Hacettepe University, Beytepe, 06532 Ankara, Turkey & 
Institute of Physics, Academy of Sciences of Azerbaijan, 370143 Baku, Azerbaijan 
e Instituto de Fisica, Universidad Nacional Autonoma de Mexico, 01000 Mexico, DF, Mexico 

Abstract 

A two-dimensional (2D) assembly of noninteracting, temperature-dependent, composite- 
boson Cooper pairs (CPs) in chemical and thermal equilibrium with unpaired fermions is 
examined in a binary boson-fermion statistical model as the superconducting singularity 
temperature is approached from above. The model is derived from first principles for the 
BCS model interfermion interaction from three extrema of the system Helmholtz free energy 
(subject to constant pairable-fermion number) with respect to: a) the pairable-fermion distri- 
bution function; b) the number of excited (bosonic) CPs, i.e., with nonzero total momenta — 
usually ignored in BCS theory — and with the appropriate (linear, as opposed to quadratic) 
dispersion relation that arises from the Fermi sea; and c) the number of CPs with zero total 
momenta. Compared with the BCS theory condensate, higher singularity temperatures for 
the Bose-Einstein condensate are obtained in the binary boson-fermion mixture model which 
are in rough agreement with empirical critical temperatures for quasi-2D superconductors. 
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1 Introduction 



Recent experiments [HJ indicate that composite bosons in ultra-cold clouds of most alkali atoms 
do indeed Bose-Einstein (BE) condense. Since Cooper pairs (CPs) of fermions (electrons or 
holes) in a many-fermion system form composite bosons in the sense of coupling to integer 
angular momentum, it is natural to consider the possible BE condensation of such pairs. The 
belief that some such condensate is central to superconductivity is more than 50 years old 
||, [| [I], ^, ||, 0]. High-T c , as well as some organic, superconductors ||] are quasi-two-dimensional 
(2D). Quasi-ID superconductors have also been found ||. BE condensation (BEC) is impossible 
in two or less space dimensions [10] for usual or "ordinary" bosons (i.e., with a quadratic energy- 
momentum, or dispersion, relation). It is however still possible to have BEC in all dimensions d > 
1 for non-interacting bosons if they obey a linear dispersion relation [11 — such as CPs moving 
in the Fermi sea. This possibility arises because the Hohenberg theorem [10], which prohibits 
BEC in 2D, relies on an /-sum rule based on the quadratic dispersion relation appropriate to 
bosons [^] moving in a vacuum. Such a linear dispersion relation for the CPs in a binary 
boson-fermion mixture model was recently found fl3| to be consistent, without any adjustable 
parameters, with the anomalous linear (quadratic) temperature-dependence above T c in the 
resistivity of optimally-doped (overdoped) cuprates whether hole- or electron-doped. For the 
observed quadratic T-dependence in overdoped samples linear-dispersion CP charge carriers are 
essential. 

Although extensive studies in the BCS-Bose "crossover" problem in superconductivity have 



spanned [14] a period of over thirty years, we note that BEC is distinct from the standard (i.e., 
zero center-of-mass momentum CPs) BCS theory condensation where only that one bosonic 
state exists. 

In this paper it is shown that in addition BEC is still possible in 2D even if the number 
of composite bosons (pairs of fermions) in a binary boson-fermion mixture is not fixed — as 
chemical/thermal equilibrium renders it coupling- and temperature-dependent — as long as the 
total number of fermions is fixed. This gives rise to an interesting statistical-mechanics problem 
irrespective of the particular mechanism for pair formation, and may have a vital application for 
superconductivity as well as for (neutral-atom) superfluidity such as in liquid 3 He p|, dilute 
mixtures of 3 He in 4 He [16|, or in trapped Fermi gases [17|. The statistical model dealt with here 
may be seamlessly linked to BCS theory, via the fermionic energy gap, when boson/unpaired- 



fermion interactions are included as, e.g., in Refs. |l^] and |T9|]. However, in these two papers 
the quadratic CP dispersion relation has been assumed. The quadratic form has recently been 
shown p0|| to apply only in the zero-density or vacuum limit when the Fermi sea disappears. 

In Section 2 we recall that at T = a 2D gas of fermions interacting via a constant pairing 
interaction in an annulus about the Fermi surface — viz., the BCS model interaction — the binding 
energy of a single pair near the Fermi surface (CP problem) decreases practically linearly with 
the center of mass momentum (CMM) of the pair for all values of the momentum below breakup, 



2 



the breakup momentum typically being only about four orders of magnitude smaller than the 
Fermi momentum. In Section 3 we discuss why the interacting many-fermion system can be 
treated as a set of independent CPs (i.e., composite bosons with fermion number two) mixed int 
with pairable fermions which are not bound into pairs, i.e., unpaired fermions. In Section 4 the 
more realistic scenario is considered of the BEC of these pairs, incorporating pair breakup beyond 
a certain CMM. Although the number of pairs is not fixed but rather strongly coupling- and 
temperature-dependent, BEC is still possible in 2D. A simple binary boson- fermion statistical 
model is introduced by constructing the Helmholtz free energy for an ideal mixture of pairable 
but unpaired fermions plus paired fermions (both zero and nonzero CMM pairs), all in chemical 
and thermal equilibrium. The latter results through extrema of the free energy in: a) the 
pairable fermion occupation probabilities; b) the excited boson numbers (nonzero CMM CPs); 
and c) the ground boson number (zero CMM pairs). In Section 5 the coupling- and temperature- 
dependence of the boson number is derived. In Section 6 the critical BEC singularity temperature 
is obtained first by ignoring the unpaired fermions in a pure boson-gas model and then exactly 
for the boson-fermion binary mixture model from a T-dependent dispersion relation derived 
and calculated numerically, and results compared with empirical data. Finally, Section 7 gives 
conclusions. 



2 Cooper-pair dispersion relation 

Consider a 2D system of N fermions of mass m confined in a square "pen" of area L 2 and 
interacting pairwise via the BCS model interaction 

Fkk , = | -V if V-P) - SwD < e fcl (= h 2 k 2 /2m), e k2 < /x(T) + hu D 
\ otherwise, 

where k = |(ki — k2) is the relative wavevector of the two particles; Vk,k' the 2D double 
Fourier integral of the underlying non-local interaction V(r, r') in the relative coordinate r = 
ri — T2] /u(T) the ideal Fermi gas chemical potential which at T = becomes the Fermi 
energy Ep = h 2 k F /2m with fop the Fermi wavenumber; Tuvd = ft 2 A;f)/2m the width of the 
annulus about the Fermi circle in which the pairing interaction is nonzero, with ujp, being the 
Debye frequency. This model interaction mimics the net effect of an attractive electron-phonon 
interaction overwhelming the repulsive interfermion Coulomb repulsions, whenever V > 0. 

If ?iK = h(ki + k2) is the center-of-mass momentum (CMM) of a pair, let Ek be its total 
energy (besides the CP rest-mass energy). The eigenvalue (CP |^I[ ) equation for a pair of 
fermions at T = immersed in a background of N — 2 inert, spectator fermions within a (sharp) 
Fermi circular perimeter of radius kp is then 

i = rv' °( kl ~ kp) 9 ^ 2 ~ kp) (2) 

V 2e k -{E K -h 2 K*/Amy 1 ' 
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where again e& = h 2 k 2 /2m, 6(x) is the Heaviside unit step function, and the prime on the 
summation sign denotes the conditions 

fci = |k + ~K| < (4 + k 2 D ) 1 / 2 and k 2 = \k-^K\<{k 2 p + k 2 D ) 1 ' 2 (3) 

ensuring that the pair of fermions above the Fermi "surface" cease interacting beyond the annulus 
of energy thickness hcop, in accordance with (|l|), thereby restricting the summation over k for 
a given fixed K. Without these restrictions @ would just be the Schrodinger equation in 
momentum space for the pair. Setting Ek = 2Ep — A^-, the pair is bound if A^ > 0, and (0) 
becomes an eigenvalue equation for the (positive) pair binding energy Ak- Our Ak and Ao 
should not be confused with the BCS energy gap A(T). 

Let A = g(Ep)V > be a dimensionless coupling constant with g(Ep) the electronic density- 
of-states (for each spin) at the Fermi surface in the normal (i.e., interactionless) state, which in 
2D is constant 



g(e) = L 2 m/2irh 2 = g. 



(4) 



The Cooper equation (Q) for the unknown quantity Ak is analyzed in Ref. [22]. For zero CMM, 
K 



0, it becomes a single elementary integral, with the familiar [21] solution 

2hu>D 



A r 



o2/X _ l 



(5) 



valid for all coupling A. For small K, it is not too difficult to extract [22] the asymptotic result 



-[1 + ^-(1 + VTT^)]hv F K + 0(K 2 ) —> A - -hv F K + 0(K 2 



F- 



(6) 



For 



where v = Op/Tp, and vf is the Fermi velocity defined through Ep = fi k F /2m 
weak coupling, A — > 0, this linear dispersion relation gives the 2D analog of the 3D result stated 



as far back as 1964 in Ref. [23], p. 33 (see also, Ref. [241 ], p. 336) but with the 2D coefficient 
2/tt of the last expression of (pi) replaced by 1/2. 



3 Justification of boson formalism 

These CP boson-like structures could be called "quasi-bosons" since their creation and anni- 
hilation operators are known not to obey the usual boson commutation relations p3|], p. 38. 
However, they do obey the Bose-Einstein distribution since the energy Ex of the CP is given 
only by the total CMM, K, but is independent of the relative momentum k. Thus, the possible 
energy states for the pair are Ek as defined in (||). The number of pairs Afc that can occupy such 
a state can take on indefinite values since there exist also indefinitely many relative momenta, 
namely 

iV K = ^Ar k)K = 0,l,2,.... (7) 

k 
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Here, A/k,K = 0, 1 is the number of pairs characterized by both k and K, and is the same number 
as that characterized by definite ki and k2, namely A/k,K = n k 1 n k 2 = 0, 1 where n^. = 0, 1 is 
the occupation number for a single fermion, these remarks all referring to singlet pairing. Much 
of all this has been known pq ] at least since 1958, albeit in somewhat different language. 

This view of an actual Cooper pair should not be confused with, say, an Anderson p3] 
phonon-like collective excitation (or modes) with weak-coupling dispersion relation — in 2D |27|] 
given by {1/ ^/2)%vfK in the long- wavelength limit, and which evolves into the plasmon when 
Coulomb repulsions between fermions are switched on. CPs here, like deuterons, carry fermion 
number two and as such are definite in number (although in the CP case this number is coupling- 
and temperature-dependent) and can thus undergo BEC. This is distinct from collective exci- 
tations which are indefinite in number. Park [28], e.g., distinguishes between "permanent" and 
"ephemeral" bosons, the latter sometimes being referred to as "quasiparticles" to distinguish 
from the former "particles" . 

For Nb ordinary bosons of mass mg and energy ek = h 2 K 2 /2mB in any positive di- 
mension, d > 0, a temperature singularity T c [29] appears in the number equation Nb = 
J2- K [e( eK ~ fMB ^ kB ' T — 1] _1 at vanishing bosonic chemical potential /i# < when the number 
of K = bosons just ceases to be negligible upon cooling. It is given by 

,2 r „_ -\2/d 



2irh z 



n B 



9d/2^) 



m B k B 

with us the boson particle density Ns/L d , and g<i/2( z ) the usual Bose integrals 

1 



(8) 



T(a) 



dx 



oo / 



z 1 e x 



1 



i=i 



(9) 



where r(<r) is the gamma function and C(c) the Riemann zeta function of order a. The last 
identification in (^) holds when a > 1 for which C(<r) < 00 > while the series 5^(1) diverges for 
a < 1, thus giving T c = for d < 2. For d = 3 one has £(3/2) ~ 2.612 so that (||) becomes 
the familiar formula T c ~ 3.31h 2 n^ 3 /mskB of "ordinary" BEC. On the other hand, for bosons 
with (positive) excitation energy ek = Ao — given approximately by the linear term in (^) 
for all K, the singularity that lead to (|l) now yields pO], for weak coupling, 



a(d)hvF 
k B 



d+l 



1/d 



(10) 



where [11] aid) = 1, 2/ir and 1/2 for d = 1, 2 and 3, respectively. Note that now T c > for 
all d > 1, which is precisely the dimensionality range of all known superconductors including 
the quasi-lD organo-metallic (Bechgaard) salts ||. This is not inconsistent with the Hohenberg 
theorem [10] that there is no broken symmetry, i.e., long-range order, in a Bose fluid for d = 
1 or 2, since this is based on an /-sum rule for bosons with a quadratic dispersion relation. 
Indeed, both (H) and (Il0|) are special cases of of the more general expression (|31fl for any space 
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dimensionality d > and any boson dispersion relation ek 
constant, given by 



C.K S with s > and C. 



Cs_ 

k B 



sT(d/2) {2ix) d n B 



s/d 



2^T{d/s)g d/s {l) ' (U) 

In what follows the number of bosons will be temperature-dependent and it is in conserving 
the fermion number that the singularity arises. As is the case for the pure boson gas, a linear 
rather than a quadratic dispersion relation will be needed to obtain BEC in 2D. This emerges in a 
statistical model of an ideal binary mixture of bosons (the CPs) and unpaired (both pairable and 
unpairable) fermions in chemical equilibrium [Q], for which thermal pair-breaking into unpaired 
pairable fermions is explicitly allowed. 



4 First-principles statistical model 

Under interaction ([!]) at any T the total number of fermions in 2D is N = L 2 kp/2iT = N± + 
N 2 and is just the number of non-interacting (i.e., unpairable) fermions N\ plus the number 
of pairable ones N 2 - The unpairable fermions obey the usual Fermi-Dirac distribution with 
fermionic chemical potential fi. On the other hand, the N 2 pairable fermions are simply those 
in the interaction shell of energy width Tuud so that 

rfl+huj D a ( e ) 

since the density of electronic states (B) is constant and the remaining integral exact. At any 
interfermionic coupling and temperature these fermions form an ideal mixture of pairable but 
unpaired fermions plus CPs that are created near the single- fermion energy /u(T), with binding 
energy A#(T) > and total energy 

E K (T) = 2fi(T)-A K (T). (13) 

This is generalizes the T = equation Ek = 2Ep — Ak introduced below (||). 

The Helmholtz free energy F = E — TS, where E is the internal energy and S the entropy, 
for this binary "composite boson/pairable-but-unpaired- fermion system" at temperatures T < T c 



is then |2| 



rfl+huJ D C \ 

F 2 = 2 deg{e)\n 2 {e)e + k B T[n 2 {e)\nn 2 {e)+{l-n 2 {e)}\n{l-n 2 {e)}]\ 



/ fl—fiUJjj 

+[2/i(T) - A (T)]N Bfi (T) 
K 



"0 f 

+ J2\[Mt)-^k(t)]n B)K (t) 

K>0 L 

+k B T [N BjK (T) In N B>K {T) - {1 + N B , K {T)} ln{l + N b ,k( T )}\ }• ( 14 ) 
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The integral term is the contribution from the unpaired fermions and runs over all levels in 
the energy shell where the BCS model interaction is nonzero, 712(e) being the average number 
of unpaired but pairable fermions with energy e; the prefactor two comes from the spin. The 
second term gives the free energy of the bosons with CMM K = since their entropy is negligible 
in the thermodynamic limit; here N B $(T) is the number of (bosonic) CPs with zero CMM at 
temperature T. The summation term represents the free energy of the bosons with nonzero 
CMM, while N ByK (T) is that with arbitrary nonzero CMM K, and the cutoff K is defined [f|] 
by Ajc = 0. The free energy F 2 is to be minimized subject to the constraint that the total 
number of pairable fermions N 2 is conserved. 

If N 2 q(T) is the number of pairable but unpaired fermions, the relevant number equation for 
the pairable (i.e., active) fermions is then 

N 2 = N 20 (T) + 2[N Bfi (T) + N Bfi<K<Ko (T)} = N 20 (T) + 2N B (T), (15) 

where Nb o<k<k (T) denotes the total number of "excited" bosonic pairs (namely with CMM 
such that < K < Ko), i.e., N b ,o<k<k (T) = 2~2o<K<K n b,k(T). Minimizing the free energy, 
subject to the constraint that (15) be a constant, is equivalent to minimizing the grand potential 



n 2 = F 2 -fx 2 N 2 . (16) 

a) Minimizing £l 2 with respect to the fermion occupation probabilities 712(e) yields 
the Fermi-Dirac distribution with fermion chemical potential fi 2 , not /j, namely 

n ^ = e^'l) + 1 ; = {ksT)' 1 . (17) 
Thus the total number of pairable (but unpaired) fermions then becomes 



N 20 (T)^2 deg(e)n 2 (e) = 2 de ^H— , (18) 



deg(e)n 2 (e) = 2 / ' 1 ' 

and should be compared with ( |l2|) for iV2 which contains only /i. Since in 2D g(e) is a constant 
(Q), (|l8|) becomes the exact expression 



N 2Q (T) = ^\n 



I _|_ g— /3(ji— (U2— Swd) 



(19) 



1 -|- e -/3{fJ.-^2+hu} D ) 

b) Minimizing ^2 with respect to the excited boson numbers Nb,k(T), K > 0, yields 
the Bose-Einstein distribution summed over all < K < Kq, namely 

K K 

N b ,o<k<k (T) = £ N b ,k(T) = E [e^ E ^~ 2 ^ - l]" 1 . (20) 

K>0 K>0 

The factor multiplying f3 in (|20|) may be rewritten as £k{T) — fi B (T), where £k(T) = Ao(T) — 
Ai<-(T) > is a (nonnegative) excitation energy as suggested by @, while [J. B (T) turns out to 
be 
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Mb(T) = 2 [ M2 (T) - fi(T)] + A (T). (21) 
This allows rewriting (|2^) in the more meaningful boson form 

Nb, 0< k<k (T) = £ [^CD-mbCT)} _ i]-i (22) 

*C>0 

where Hb{T) is clearly the bosonic chemical potential associated with the entire binary mixture. 

c) Finally, minimizing 0,2 with respect to the number of zero CMM (or, "ground 
state") bosons Nb,q(T) gives 

2MT) - n{T)\ + A (T) =0 (0 < T < T c ), (23) 

valid only in the stated temperature range as Nb,o(T) is negligible for all T > T c . However, 
in view of ( ^T|) this implies that Hb(T) = for all < T < T c — which is precisely the BEC 
condition for a pure boson gas, even though one now deals with a binary boson-fermion mixture. 



5 Boson number 



To determine Nb{T) from (15) we need (^) which with (^) reduces to 



N 20 (T) = ^ln 



1 + e -/3{A (T)/2-^ I5 } 
1 _|_ e -/3{A (T)/2+^ I5 } 



(0 < T < T c 



(24) 



At T = two distinct coupling regimes emerge by inspecting (|24|): a) for Ao/2 < Hujd or, from 
© for A < 2/ In 2 ~ 2.89, we have that iV 2 o(0) = 2g(^)(fiuj D - A /2); while b) for A /2 > %u D 
(or A > 2.89) ^20(0) is identically zero. Hence, the number of bosons Nb(0) at T = from ( |T5| ) 
is just Nb(0) = o[-^2 — -^20(0)]. Using (12) for N2 the fractional number of pairable fermions 
that are actually paired at T = 0, namely 2Nb(0)/N2 = 1 — A^oW/A^, becomes simply 



A /2hio D = (e 2 / A - I)" 1 — ► e^ 2 / A (for A < 2/ In 2 ~ 2.89) 
2N B (0)/N 2 = { ^0 ~ (25) 

1 (for A > 2/ In 2 ~ 2.89). 

This fraction is plotted against coupling A in Fig. 1. Since N B {0) = ^A for A < 2.89, only 
those fermions in an energy shell of width Ao/2 around the Fermi surface actually pair at T = 0, 
while for A > 2.89 all pairable fermions actually pair up since then N B (0) = ghu D = \N 2 . This 



result contrasts sharply with the "heuristic model" |3l|, Eq. (22), where 2Nb(0)/N2 = 1 for all 
coupling, and is more in line with BCS theory which implies, in any d, a coupling-dependent 
fraction estimated (Ref. § p. 128; see also |1) to be [g{E F )2A/2g(E F )huj D } 2 = (&/Tiuj d ) 2 = 



(sinhl/A) 2 — * 4e 2 / A , where A = fiujo/ sinh(l/A) (again, not to be confused with the CP 

A — >0 
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binding energy Ao) is the T = BCS energy gap for the same BCS model interaction (||) used in 
this paper; this is graphed as the long-dashed curve in Fig. 1 and is seen to be much larger than 
(|25| ) for fixed A. The breakdown of BCS theory for BCS model interaction couplings larger than 
A ~ 1.13 is clear both because: a) the alluded fraction cannot exceed unity; and b) physically, 
if the fermionic energy gap A > Tuod no pairable fermions are available at all. This breakdown 
is indicated by the short-dashed curve in Fig. 1. (A strong-coupling many-body model differing 
from that of BCS theory but based on the BCS model interaction has been solved by Thouless 




Also displayed in Fig. 1 are two finite-temperature results for 2Nb{T)/N2 = 1 — N2$(T) / N2 
which are obtainable from fl24j ) for any T provided one knows Aq(T) for any T > 0. For T > 0, 
the 8(ki — kp) = 9(ek 1 — Ep) in (Q) becomes 1 — n(^ 1 ), where = (e^ fc i + 1) _1 is the 

Fermi-Dirac distribution with ^ = e/ Cl — fJ,(T), with the ideal fermion gas chemical potential 
jU(T) in 2D being given exactly by 

f,(T) = /3- 1 He PEF -l)^ Ep. (26) 

Note that /jl(T) decreases monotonically with temperature from its maximum value of Ep but 
does not turn negative until T = Tpj In 2 ~ lAATp so that the BCS model interaction 
which requires fi(T) to be nonnegative, will not break down (i.e., become meaningless) over the 
entire range of temperatures relevant in this paper, see Fig. 4 below. Similar arguments hold 
for 0(&2 — kp). Since k\ = k2 implies that £fc a = £& 2 , (Q) then leads to a simple generalization to 
finite-temperature of the K = CP equation, namely 

1 = A / de(e~« + iy 2 [2t + Ao(T)]- 1 . (27) 
Jo 

Its numerical solution for Aq(T) is illustrated in Fig. 2 for A = gV = 1/2 and v = Tilod/Ef = 
0.05. Note that if one assumes a T* such that Ao(T*) = 0, the resulting integral in ( |27| ) diverges 
and the equation can only be satisfied for A = 0; thus, there is no temperature T* at which 
"depairing" will occur for any fixed v and any nonzero A. 



6 Critical temperature 

Neglecting the background unpaired fermions and modeling our system as a pure boson gas 
of CPs but with temperature-dependent number density n_e(T), one converts the explicit T c - 
formula ([h]) into an implicit one by allowing to be T-dependent. For d = 2 (pi) becomes, 
since 8,(1) = C(2) = ir 2 /6, 

4^ Hvp 



T c = ^^n B (T c ). (28) 

This requires ub (T) = N B (T)/L 2 which in turn requires (J24), along with Aq(T) as determined 
from @, and is given by the expression 2N B (T)/N 2 = 1 - N 2 o{T)/N 2 . Solving (||) self- 
consistently with A = 1/2 gives the remarkably constant value T c /Tp ~ 0.004 over the entire 
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range of v = Tllod/ Ep values 0.03 — 0.07 typical [|35| of cuprate superconductors. On the other 
hand, the BCS formula Tf cs ~ 1.139 £) e- 1 / A with A = 1/2 gives T c /7> = 0.005, 0.008 and 
0.011 for v = 0.03, 0.05 and 0.07, respectively. Clearly, both sets of predictions are somewhat 
small compared with empirical cuprate values of T c /Tp that range |3(| from 0.01 — 0.1. 

To obtain the exact critical temperature without neglecting the background unpaired fermions, 
one needs the exact CP excitation energy dispersion relation £k(T) = Aq(T) — Ak(T) which 
is neither exactly linear in K nor independent of T. To determine A_r-(T) we need a working 
equation that generalizes Ref. |2^] for T > via the new CP eigenvalue equation ( ]2^ ) . Because 
of symmetry, see Fig. 3, one can restrict the angle 9 to the interval (0,tt/2) where k\ > k2, i.e., 



to quadrant I. Recalling (13), in d-dimensions (J2j) becomes 



i=v(—Y r'dk [i-^foj] [!-»(&)] m 

\2ttJ J h 2 (k 2 -k 2 )/m + A K (T) + h 2 K 2 /4m' V ' 

Here k^ is such that /u = h~ 2 k 2 l /2m and becomes kjr as T — > 0, while fco is such that ?k<j£) = 
h 2 k 2 ~ ) /2m. The prime on the integral sign now denotes the restrictions 

4 = |k-±K| 2 = k 2 -kK cos 6 + \K 2 > A; 2 , (30) 
fc 2 = |k+±K| 2 = k 2 + kKcos6 + \K 2 < k 2 +k 2 D . (31) 

In Fig. 3 the darkest shading corresponds to these (BCS model interaction) restrictions. The 



conditions ( p0[) and (^Ij) can be studied separately but must be satisfied simultaneously. If 



K < 2y 'fc2 _ fc2 )) (|oD an d (|i|) are equivalent to 

(A; 2 - \K 2 sin 2 0) 1 / 2 + ±K cos 9 < k < [(k 2 + k 2 D ) - \K 2 sin 2 9} 1 ' 2 - \K cos 9 . (32) 



Note that for K > ^ k"^ + k 2 D — k 2 , — k 2 D there exists a minimum value # m j n of given by 

K^2{2kl + kl)-K 2 



k 2 

cos (9 min = g _ , (33) 



while m j n = for i^T < + 4 — \Au — ^i>- We introduce the dimensionless variables 



2(fc 2 +4)1/2 - 1 ' S-^' ^- jBf ' V -T F ~k\ 



-))'• <vz<' J-t k 2 

with ks@D = ftc^D = h 2 k 2 - ) /2m and UbTf = Ep, where &_b is Boltzmann's constant. Recall 
the d = 2 constant expression (Q) for (7(e), the restrictions (^), and that for K > and T > 
the step functions in @ 0(A;i )2 - fcj?) = 6>(|±K ± k| - fcj?) become [exp{-/3[?i 2 (±K ± k) 2 /2m - 
jU(T)]} + 1] _1 — but with 2efc in (||) replaced by ek 1 + £fc 2 , -Ef by n(T) and by Ak(T). One 
finally arrives at a working equation for the binding energy A#-(T) that generalizes Eq. (18) of 
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Ref. j22j, namely 



* h mm Juin(e) 2£ 2 + 2(1 + v)k? — 2 + A K (T) 

x[l + exp{-/3[£ 2 + (l + z/)K 2 -2 v / m ; <cos6<- l]}]" 1 , (35) 



where i/ = huj D /fi, £, min (9) = ^fl + v kcos0 + J I - (1 + ^)k 2 sin 2 6*, £ TOaa; (0) 



VI + ^ kcos6> + x /(1 + z/)(1-k 2 sin 2 6») and 



ff f if 2 K <l-y/(l-u)/(l+u), 

mm \ cos _1 (z^/{4Vl + v K^J\ + u/2 - (1 +v)k 2 }) otherwise. 
In (j35|) we have introduced the more general dimensionless quantities £ = k/k^, A K (T) = 



A K (T)/fi, where T = /cbT/^ or $ = fj,/3, and k = K/2^Jkl + fcf,. 

To obtain the critical temperature from the finite-temperature dispersion relation, besides 
solving (f2§) for A K (T), one needs (0), @ and At T = T c both N B ,o(T c ) ~ and 

fJ-B(T c ) — so that one gets the implicit T c -equation for the binary mixture gas 



l = ^ln 

v 



1 + e -{Ao(T c )/2-»/}/f e 



8(l + i/) r K o(Tc) K 

+ V 7 K e [A (f c )-A K (f e )]/f c _ 1 - ! " >(>! 



1 + e -{A (T c )/2+y}/T c 

This must be solved numerically for the exact T c for each A and v in conjunction with (p7|) for 
Ao(T) and ( |35| ) for both A K (T) and Ko(T c ). Results for A = 1/2 are shown in Table 1 and 
Fig. 4 for a range of v values typical |35[] of cuprates. We have taken T^/Tp ~ 1, a very good 
approximation up to the highest temperatures dealt with. For example, from Fig. 4 the highest 
T c /Tp ~ 0.14 already gives T^/Tp ~ 0.9999 from (p6|), while for smaller T c /Tp the values of 
T^/Tp are even closer to 1. The T c resulting from the exact dispersion relation for T = 
(dot-dashed curve) is somewhat higher than the exact result (full curve) but lower than that 
using the linear approximation for Ak (T) (dotted curve) . It is also clear that the effect of using 
the exact or linear (in K) cases dominates the effect of the dispersion relation T-dependence. 



For cuprates d ~ 2.03 has been suggested [37] to be more realistic as it reflects inter-CuO-layer 
couplings but our results in that case would be very similar to those reported here for d = 2. 

Thus, for v = 0.05 the exact T c is seen to be about 46% lower than the heuristic result 
found in Ref. [3l[] , Eqs. (|l^) and (p3|) . It is curious that all results depend very weakly on the 
T-dependence of the CP binding energy Ak{T), in spite of its being substantial throughout the 
temperatures spanned in this paper, as seen in Fig. 2. 

We defer study of the condensate fraction Nb,q(T) /Nb(T) below T c and merely surmise that 
it may ultimately help explain the apparent absence [ |38| , |39"[ in cuprates of the Hebel-Slichter 
peak of nuclear-spin (NMR) relaxation rates vs temperature for < T < T c . Such a peak, 
originally seen M3] in aluminum, is perhaps the most stringent and qualitatively convincing 



experimental test of BCS theory (Refs. [23], p. 71 and [41|, p. 79 ff). Besides cuprates, it 
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is also absent [42] in several quasi-lD Bechgaard ||] and in several quasi-2D (ET) organic salt 
superconductors . 



7 Conclusions 

A simple statistical model treating CPs as non-interacting bosons in thermal and chemical equi- 
librium with unpaired fermions is proposed. The model gives rise to a boson number that 
is strongly coupling- and temperature-dependent. Since the CP dispersion relation is approx- 
imately linear, it exhibits a Bose-Einstein condensation of zero-CMM pairs at precisely two 
dimensions. Exact transition temperatures based upon the exact CP dispersion relation are in 
reasonable agreement with empirical cuprate data. 

Needless to say, further corrections are yet to be included in the present simple binary 
mixture boson- fermion model, e.g., i) realistic Fermi surfaces, ii) Van Hove singularities [O] 
or other means of accounting for periodic-crystalline effects, as well as iii) the all important 
d-wave interfermionic interaction, iv) the boson-fermion interaction and v) residual interbosonic 
interactions. As to the latter, also generally neglected in BCS theory, if the lowering [44] of T c 



in liquid 4 He by about 29% with respect to the ideal Bose gas BEC T c is any guide, interbosonic 
interactions will also lower T c in a more realistic picture. As to the boson-fermion interaction, it 



is precisely this ingredient that enabled T.D. Lee and coworkers [18], and Tolmachev [19] more 
generally, to link BCS and BEC through a relation stating that the BE condensate fraction is 
proportional to the (BCS-like) fermionic gap A(T) squared. 
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Table 1: Critical temperatures T c /T^ for A = 1/2 depicted in Fig. 4 according to (13(f). The 
exact result is compared with the linear- in- K approximation for both Ak(T) and A^(0) in 
order to test sensitivity of a temperature-dependence of the CP binding energy for nonzero K. 



V 


linear approx. with A^(T) 


linear approx. with A^-(O) 


Exact 


0.03 


0.078 


0.068 


0.065 


0.04 


0.089 


0.079 


0.075 


0.05 


0.100 


0.088 


0.084 


0.06 


0.109 


0.096 


0.091 


0.07 


0.117 


0.104 


0.098 



Figure Captions 

Figure 1. Fractional number of pairable fermions that are actually paired, at three different 
temperatures, vs. coupling A for the present first-principles model (p5|) and estimated for BCS 
theory at T = as explained below (^). The number of pairable fermions with the BCS model 
interaction used is just (12); all of them are actually paired at T = in the heuristic BEC 
model, Ref. @ Ec l- ( 23 )- 

Figure 2. Temperature dependence of K = CP binding energy Ao(T) obtained numeri- 
cally from ( p7| ) for A = 1/2 and v = 0.05. Note that when T = oo ( p7| ) is analytical for Ao(oo); 
the latter then turns out to be about 10 -8 , so that the curve saturates from above to this value 
at T = oo. 

Figure 3. Cross-section of overlap "volume" in momentum space (darkest shading) where 
the tip of the relative wavevector k (for two fermions with wavevectors ki and k2) must point 
for the attractive BCS model interaction (||) between them to be nonzero and form a Cooper 
pair of CMM magnitude TiK. 

Figure 4. Critical BEC temperature T c in units of Tp, resulting from (|36|) for A = 1/2 for 
varying v = Hujd/^ — @d/Tf: with no approximations (full curve); using Ak(T) evaluated 
at T = (dot-dashed); using the linear- in- K approximation for A^(T) (dotted). The dashed 
straight line is the BCS formula T c ~ l.l?>Q D e l/2X for A = 1/2. The very lowest full horizontal 
line is the solution of the implicit T c -equation (|28|) for the pure boson gas for v = 0.03, 0.05 and 
0.07. Cuprate data are taken from Ref. |3£ 
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